Role of confined phonons in thin film superconductivity 



E. H. Hwang and S. Das Sarma 
Center for Superconductivity Research, 
Department of Physics, University of Maryland, College Park, Maryland 20742-4111 

M. A. Stroscio 
U. S. Army Research Office, P.O. Box 12211 
Research Triangle Park, NC 27709-2211 
(February 7, 2008) 

We calculate the critical temperature T c and the superconducting energy gaps A n of a thin film 
superconductor system, where A n is the superconducting energy gap of the n-th subband. Since the 
quantization of both the electron energy and phonon spectrum arises due to dimensional confine- 
ment in one direction, the effective electron-electron interaction mediated by the quantized confined 
phonons is different from that mediated by the bulk phonon, leading to the modification of T c in 
the thin film system. We investigate the dependence of T c and A„ on the film thickness d with this 
modified interaction. 
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Superconductivity in thin films has been studied for 
the last four decades. The phenomenon of thin film su- 
perconductivity has its p<pwn specific peculiar features. 
In early investigations^ the effect of the film size on 
the superconducting transition temperature, T c , for thin 
films was investigated. Experiments have shown a mono- 
tonic increase of the critical temperature, T c , with de- 
creasing film thicknessa. From the theoretical point of 
viewcTEl, the shape resonances and the strong thickness 
dependence of T c are the characteristic features of thin 
film superconductivity. The size quantization of the 
transverse motion of the electron in the film leads to 
an increase of T c with decreasing film thickness, arising 
essentially from an enhanced effective BCS pairing in- 
teraction. The resonance effects are manifest each time 
one of the n two-dimensional (2D) subband energy levels, 
E n (d), of a film with thickness d, (for 'transverse' motion 
perpendicular to the film faces or along the confinement 
direction), passes through the Fermi surface as the thick- 
ness d is varied. In previous calculationsau, the phonon 
modes were assumed to be the same as in the bulk mate- 
rial and only the one-dimensional quantum confinement 
effects of electrons were considered; i.e., the superconduc- 
tivity in the thin film was considered to be arising from 
the attractive electron-electron interaction mediated by 
bulk phonons interacting with subband-quantized elec- 
trons. However, phonons in thin films ('slab' phonon 
modes) also have specific characteristic features by virtue 
of definite boundary conditions imposed by confinement 
in thin films. The phonon dispersion in a thin film under- 
goes substantial modification compared with the bulk, 
and a quantization of the phonon spectrum occursErQ. 
The quantization of the phonon spectra has many ef- 
fects on thin film electronic properties, which have been 
extensively studied, particularly in the context of semi- 



conductor quantum well structures. It is the puxpose of 
this paper to reconsider the previous calculations!!! of thin 
film superconductivity taking into account confined size- 
quantized phonons in the thin film. We find the resonant 
shape of the superconducting transition temperature, T c , 
arises from both the quantum electronic confinement and 
phonon confinement. _ 

Following the BCS theorya the critical temperature 
T c ~ 1.14o;£)exp(— 1/iVbVb) for a bulk BCS supercon- 
ductor, where wp> is the phonon Debye frequency, Nq the 
electronic density of states at the Fermi energy, and Vb 
the effective attractive electron-electron interaction me- 
diated by phonon exchange. An increase of NqVq im- 
plies an increase in T c . For a thin film with thickness 
d, the density of states at the Fermi energy develops 
quantum structure due to the confined 2D subbands: 
Nfum = (2u - l)m/47rd = N 3D (2v + l)-n/2dk F , where 
v = 1, 2, ... is the occupied subband index. Thus, when 
we use the bulk electron-phonon coupling constant J 
unmodified by any phonon size confinement corrections, 
which is identical for all subbands, the critical temper- 
ature for thin films depends on the film thickness, and 
decreases exponentially as the thickness increases. How- 
ever, for a fixed bulk electron density (or, hp), as the 
film thickness increases the higher quantized subbands 
are occupied by Cooper pairs, so that the critical tem- 
perature of the film jumps to a higher value due to the 
higher v, arising entirely from the jumps in the density of 
states as the effective Fermi level moves through higher 
values of the subband index v. This implies that with 
increasing thickness of the film the critical temperature 
of the film exhibits resonance features. In real thin films, 
therefore, the electron-phonon coupling constant is differ- 
ent from that in the bulk material due to the quantiza- 
tion of the phonon dispersions, and therefore the simple 
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resonance scenario (discussed above) must be modified. 
When we consider the confined phonons in the thin film, 
the electron-phonon coupling constant is different for dif- 
ferent subbands, and the effective coupling strength de- 
creases with increasing subband index; that is, higher 
subbands have progressively weaker electron-phonon cou- 
pling constants. This change of the coupling constant 
gives rise to features in the superconducting energy gap 
of thin films, which have not earlier been considered in 
the literature. In this paper we calculate the critical tem- 
perature, T c , and superconducting energy gaps, A„, with 
the modified electron-electron interaction mediated by 
the confined slab phonons. The results for the calcu- 
lated critical temperature obviously depend on the slab 
phonon dispersion, which in turn depends on the bound- 
ary conditions used in the phonon confinement model. 
We have followed here the simple approach of Thomp- 
son and Blattld, where the boundary has been treated by 
an infinite wall, and we use wavefunctions which vanish 
at the boundary. Although this is a highly simplified 
phonon confinement model it has the virtue of being an- 
alytically tractable — one could systematically improve 
upon this model using our theory as the starting point. 
We expect our results to be well-valid qualitatively. 

We assume that our superconducting film with a fi- 
nite width, d, is confined in the z direction by an infinite 
square-well confinement potential applied at z = and 
z = d. We choose the same infinite (one dimensional) 
square- well confinement for both electrons and phonons, 
as would be appropriate for a free-standing thin film. 
Using periodic boundary conditions in the x and y di- 
rections with periodicity distance L, we have the one- 
electron wave function and spectrum 



0k,n(r, z) = u n (z)exp(ik • r)/L, 
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Here, k = (k x ,k y ), r = (x,y) are the 2D wave vector 
and position vector in the plane of free motion, and by 
solving Schrodinger equation for z-direction with confine- 
ment potential we have 



u n {z) = (2/d) 1 / 2 sm(k n z), 



E n = 
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where k n — nn/d with n = integer. Thus, confinement 
in the z direction leads to the quantization of electron 
energy levels into different subbands. In addition to the 
quantization of electron energy, we take into account the 
modification in the thin film phonon dispersion -arising 
from the quantization of the phonon spectrum&Ll. The 
quantization of the phonon spectra leads to the change of 
the conventional electron-phonon interaction. The spe- 
cific expression for the electron-phonon coupling in the 
thin film can be obtained on the basis of the general de-i 
formation potential electron-phonon interaction theoryt2l 



where 



H cp = D I d 3 r$> 1 (r)Vw(r)^(r), 



*(r) = -j= J2 u n{z)e lk ' r c kr , 
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A is the 2D the destruction operator for an 

electron with momentum k in the n-th subband, w(r) 
is the lattice displacement vector, D is the deformation 
potential coupling constant (which we uncritically as- 
sume to be the bulk value since not much experimen- 
tal information is available). The displacement vector 
w(r) can be obtained by solving equations of motions 
of the elastic continuum with the appropriate boundary 
conditiona!ffl. We skip the details of this elastic contin- 
uum theory (which involves substantial algebra) for the 
sake of brevity 

Within the framework of the BCS theory, Cooper 
pairs are produced consisting of electrons with oppo- 
site k and identical n in a thin superconducting film, 
i.e., the electron pairs are produced between (k, n, a) 
and (— k, n, — a). The direct effective electron-electron 
interaction due to the exchange of virtual confined slab 
phonons becomes 



V nm (q) 



^Mr(q)Di(q)M[' 
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where M™ m is the electron-confined phonon matrix el- 
ement and Di the phonon propagatora, and the I sums 
cover all phonons with energy less than the cutoff en- 
ergy, cod- The maximum value of I contributing to the 
sum over the slab phonon mode I in Eq. (|^), l ma x, is 
given by the condition l ma x — {d/Tr)(huj]j/c), where c is 
the velocity of the phonons. (We assume that all slab 
phonons have the same velocity because information on 
the slab mode velocity is not experimentally available.) 
We then get, 
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where the confinement form factor is given by 
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Eq. (^) has the following selection rule: for |n±n'| = even 
(odd) only odd (even) Z's are allowed in the sum over slab 
modes I. The slab phonons with odd / are symmetric and 
the ones with even / are antisymmetric with respect to 
reflection through z = d/2, i.e., under the transformation 
\z — d/2\ — > \d/2 — z\. Since the electron wave functions 
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are either symmetric or antisymmetric, the couplings be- 
tween two subbands of the same symmetry involve sym- 
metric quantized phonons, while couplings between two 
subbands of different symmetry involve antisymmetric 
phonons. (Note that this simplicity will be lost if parity is 
not a good quantum number as it is in our simple infinite 
square well model, but would not be under an asymmetric 
confinement.) We see that V nm decreases with increas- 
ing /, since the transition of electrons to higher subbands 
cannot be induced by phonons with small momenta q. 
Thus, the components V nm of the interaction matrix form 
a monotonically decreasing sequence with increasing sub- 
band index. When all the confined phonons contribute 
to V nm (i-e., I max —> oo) we recover the bulk phonon me- 
diated resultsEl, i.e., V nm = — (J/Ad)(l + \8 nn '), as we 
should. As a result of the confined phonons, a supercon- 
ducting condensate of Cooper pairs can be produced in 
a given subband, first because of the attraction due to 
the electrons in the same subband and second because 
of transitions from other subbands contributing to the 
condensate. 

In the presence of a number of subbands the reduced 
BCS Hamiltonian of the system is given by 



k,n.cr 
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non-linear coupled subband matrix equation. The critical 
temperature is given by 



T c = I.Ulod exp 
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where X\ n — A n /Ai, the ratio of the n-th subband en- 
ergy gap to the lowest subband energy gap. For any 
given finite width, d, of the slab, only a finite num- 
ber of eigenvalues, £ n , contribute; values of £„ in ex- 
cess of fi + hup make vanishing contributions, because 
then all the £k's lie outside the interaction region. Thus, 
the summation in Eqs. ([l0|) and ([ll]) is only over all 
the occupied subbands. For a fixed electron density we 
can find the maximum value of the n (or, the highest 
occupied subband) from the chemical potential. The 
number density no = N/V is related to the chemi- 
cal potential by the relation, N — nit n , where 

Tik., n = [exp(— £ n (k)/kBT) + is the Fermi distribu- 
tion function. At T = OK we have 

H = (irdh 2 /isM) Lo + + \){» + 1) j , (12) 

where v is the maximum value of occupied subband n 
and is given by the integral value of the expression, 
v = dkp/n. If we write Eq. ( pT| ) in the usual BCS 
form, T c = 1.14wd exp(— 1/Nfu m J e ff), the effective in- 
teraction parameter can be written as 



where 



'kner 



is electron creation operator in the n-th sub- 



band with spin a, £ n (k) = e„(k) — /i the electron energy 
in the n-th subband measured from the chemical poten- 
tial /j,, and Vnm the attractive interaction between n-th 
subband and m-th subband mediated by the confined 
phonons in the film. 

In the BCS theory (which is what we utilize) the gap 
function has the same energy cutoff Uu>d as the interac- 
tion. In the weak-coupling approximation we have the 
the superconducting energy gap for n-th subband given 
by the gap equation 
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where E n — + A^) 1 / 2 . The nondiagonal terms in 
the sum reflect the possibility of the transition of the 
electron pair from one subband into another as a result 
of interaction with confined phonons. Integration over k 
gives the gap function of the subband n at T = K 
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where a nn ' — V nn ' / J ■ If all confined phonons contribute 
equally to the electron-electron interaction, i.e., (cx nm — 
1 + S nn i/2), then we recover the results of ref. [j]. With 
the coupled subband interactions Eq. (10) becomes a 
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As d — * oo, J e ff — * J because a\ n {d — > oo) = 1 + 5\ n /2 
and X\n(d oo) = 1, and Nfu m —* N^d- Thus, we 
recover the three dimensional results as we expect. 

In Fig. 1 we show the calculated superconducting en- 
ergy gaps of each subband and the critical temperature 
as a function of the film thickness d up to the 4-th sub- 
band occupation. In this figure we use the following pa- 
rameters: Debye energy Tiujd — WOK, electron density 
n = 2 x 10 22 cm" 3 , and N 3D J = mk F J/{2n 2 ) = 0.3. 
In previous calculationsBH the energy gaps of different 
subbands were find to be the same by virtue of the bulk 
phonon approximation—The shape resonance feature in 
the earlier calculationsoo arise only from the effective 
2D density of state of the film as the chemical potential 
passes through different subbands. In our calculation 
we find that the energy gaps are different for different 
subbands since the effective coupling strength depends 
explicitly on the occupied subbands. In addition, the 
resonance structures in our results (the sharp maxima in 
Fig. 1) arise from both the thickness dependent density 
of states Nfum and the effective interaction parameter 
J e ff- The energy gap is maximum for the lowest sub- 
band and decreases in the higher subbands for a fixed 
film thickness. In Fig. 1(b) we also compare our results 
with Ref. ^ for equivalent parameter values. Our results 
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exhibit more resonance features and has much lower T c 
(typically, about half around the maxima) than that of 
Ref. 0. Note that any inhomogeneity on the microscopic 
scale (quite unavoidable in real thin films of 5-15 A thick- 
ness) will considerably suppress the resonance features of 
Fig. 1, and any enhancement in A n or T c may remain 
unobservable unless the films are microscopically of uni- 
form thickness. The reduction of the critical temperature 
in our calculation can be explained by the enhancement 
of the effective interaction parameter in the slab phonon 
model. The present work can easily be extended to in- 
clude more realistic boundary conditions! 2 ], but one then 
needs to resort to numerical work right from the begin- 
ning, losing much of the essential qualitative physics of 
the phenomenon. We believe that the basic physics dis- 
cussed in this paper and the qualitative features of our 
results shown in Fig. 1 transcend our specific model, and 
should be valid in any BCS type superconductivity in 
thin films. 
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FIG. 1. (a) The calculated superconducting energy gaps, 
A n , of each occupied subband n, and (b) the critical temper- 
ature for an electron density n — 2 x 10 22 cm~ 3 as a function 
of the film thickness d. In (a) the highest curve is the energy 
gap of the lowest subband (n = 1), and the second highest 
curve is that of the first excited subband (n = 2), and so 
on, with four occupied subbands. In (b) the solid line in- 
dicates our result with the interaction mediated by confined 
slab phonons and the dashed line from ref. [1] corresponds 
to the bulk phonon result. The arrow in each figure indicates 
the purely bulk result. 
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